Spin-stiffness of anisotropic Heisenberg model on square lattice and possible 
mechanism for pinning of the electronic liquid crystal direction in YBCO 
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Using series expansions and spin- wave theory we calculate the spin-stiffness anisotropy p 3X /p ay 
in Heisenberg models on the square lattice with anisotropic couplings J x ,J y . We find that for the 
weakly anisotropic spin-half model (J x ss J y ), p sx /p S y deviates substantially from the naive estimate 
p sx / Pay ~ Jx/Jy We argue that this deviation can be responsible for pinning the electronic liquid 
crystal direction, a novel effect recently discovered in YBCO. For completeness, we also study the 
spin-stiffness for arbitrary anisotropy J x /J y for spin-half and spin-one models. In the limit of 
Jy/Jx — * 0, when the model reduces to weakly coupled chains, the two show dramatically different 
behavior. In the spin-one model, the stiffness along the chains goes to zero, implying the onset 
of Haldane-gap phase, whereas for spin-half the stiffness along the chains increases monotonically 
from a value of O.I8J2; for J y / J x = 1 towards 0.25Ja; for J y /J x — -> 0. Spin-wave theory is extremely 
accurate for spin-one but breaks down for spin-half presumably due to the onset of topological terms. 

PACS numbers: 



This work is motivated by the recent discovery-^ of 
the electronic liquid crystal in underdoped cuprate super- 
conductor YBa2Cu3 06.45- The electronic liquid crystal 
manifests itself in a strong anisotropy of the low energy 
inelastic neutron scattering. The liquid crystal picture 
implies a spontaneous violation of the directional sym- 
metry: the "crystal" can be oriented either along the 
(1,0) or along the (0,1) axes of the square lattice. The 
YBa2Cu3C>6.45 compound has a tetragonal lattice with 
tiny in-plane lattice anisotropy, a* /b* ~ 0.99. This tiny 
anisotropy is sufficient to pin the orientation of the elec- 
tronic liquid crystal along the a*-axis. As a result the 
low energy neutron scattering^ demonstrates a quasi- 
1D structure along a*. 

To understand the pinning mechanism of the electronic 
crystal, in the present work, we study the anisotropic 
Heisenberg model. We calculate the in-plane anisotropy 
of the spin-stiffness and demonstrate that this is strongly 
enhanced by quantum fluctuations. We argue that the 
enhancement is sufficient to provide a pinning mechanism 
for the initially spontaneous orientation of the electronic 
liquid crystal and suggest a specific mechanism for the 
pinning. 

The anisotropic Heisenberg model has previously at- 
tracted a lot of theoretical interest&^&i. However, most 
theoretical studies have focussed on the regime of strong 
anisotropy, where the system reduces to one of weakly 
coupled spin-chains, and the most significant issue there 
is that of dimensional crossover and the onset of long 
range antiferromagnetic order. To the best of our knowl- 
edge the anisotropy of spin-stiffness has not been studied 
before. This is an important theoretical problem in it- 
self and therefore we extend our study to the case of 



arbitrary strong anisotropy. We consider both spin-half 
model, where in the limit of strong anisotropy we come to 
the situation of weakly coupled Heisenberg S=l/2 chains, 
and also spin-one model, where in the limit of strong 
anisotropy we come to the situation of weakly coupled 
Haldanc chains. 



Our series expansion results show that the spin- 
stiffness indeed behaves very differently in the two cases. 
For spin-one, the stiffness along the chains vanishes at 
an anisotropy ratio of J y /J x ~ 0.01. Self-consistent spin- 
wave theory remains highly accurate in this case all the 
way down to the transition. On the other hand, for spin- 
half, series expansions show that the stiffness along the 
chains increases from 0.18J X in the isotropic limit to- 
wards the known£ ID result of 0.25^ as J y — > 0. In 
this case spin-wave theory clearly breaks down with in- 
creasing anisotropy, presumably due to the onset of Berry 
phase interference^. 



The structure of the paper is as follows: in Scction|T]we 
calculate the spin-stiffness using series expansions. This 
is probably the most accurate method that is valid from 
small to very large anisotropy. In Section |TT] we calcu- 
late the same spin-stiffness using spin-wave theory. This 
method is valid as long as one is not close to ID limit. 
In Section HTT1 we discuss the application of our results to 
the explanation of the electronic liquid crystal pinning 
in YBa2Cu30g.45. Finally, in Section llVl we draw our 
conclusions. 
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FIG. 1: (Color online) Series Expansion (SE) and Spin-Wave 
(SW) spin-stiffness of the spin-half Heisenberg model on the 
anisotropic square lattice along the x and y axis as a function 
of the anisotropy J y /J x . The dotted line shows the value of 
the spin-stiffness p sx in the ID limit of the model (J y /J x = 0). 
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FIG. 2: (Color online) Series Expansion (SE) and Spin- Wave 
(SW) spin-stiffness of the spin-one Heisenberg model on the 
anisotropic square lattice along the x and y axis as a function 
of the anisotropy J y /J x . In the region 0.01 < J y < 0.1 the 
data points for the SE p ax have been fitted to the curve p sx — 
( J y — 0.01)", with v = 0.7±0.1. The dotted line represents the 
spin-stiffness in the 2D isotropic limit of the model (J y /J x — 
!)• 
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FIG. 3: (Color online) Series Expansion (SE) and Spin- 
Wave (SW) spin-stiffness anisotropy p ax /p ay as a function 
of J x /J y — 1 for the spin-half Heisenberg model on the 
anisotropic square lattice. The data points have been fitted 
to the linear function given in Eq. (J5J) yielding use — 1.8 and 
Ksw = 1-3 



the ground state energy of the system under an applied 
twist along one of the axe s 10 ! 11 . In general, it can be 
decomposed into a sum of two parts, a paramagnetic part 
and a diamagnetic part. For the anisotropic model, one 
can define two different twists p sx and p sy depending on 
whether the twist is applied along the x or the y axis. 
Following Ref. llolflTl . the diamagnetic component of the 
twist for p sy is given by the expression 



dia 
Hsy 



-J <r ^ z S" 2 

J V ^ °r J r+y 



S'pSp +y > 



(2) 



where angular brackets denote expectation value in the 
ground state of the Hamiltonian in Eq. {l}. The param- 
agnetic term is given by the equation 



jpara op 

Psy - ltj S , 



(3) 



where Eg is the coefficient of the 9 2 term in the ground 
state energy per site of the Hamiltonian in Eq. |T]) with 
a perturbation 



n/para t n \ ^ qx ( qz qz \ 

n ~ J y° 2_*i r y°r+y D r-y) 



(4) 



I. HAMILTONIAN AND SERIES 
CALCULATION 

We consider Antiferromagnetic Heisenberg model on a 
square-lattice, with spatially anisotropic exchange cou- 
plings given by the Hamiltonian 

W = J X ^ ' Sf ■ Sf+& + Jy ^ ' Sf> ■ Sf-\-fi , (1) 

r r 

where the sum over r runs over all sites of the square- 
lattice. Spin-stiffness can be defined by the change in 



In order to calculate these quantities we introduce an 
Ising anisotropy^ by scaling all XY parts of the exchange 
interactions by a factor A. Then p sx and p sy can be 
calculated as a power series in A for any value of the 
coupling anisotropy. Scries expansions for selected values 
of the anisotropy for the spin-half and spin-one model are 
given in Tableland Table [Hi respectively. 

The series are analyzed by Integrated Differential Ap- 
proximants (IDA)^. Before the analysis, a change of 
variable of the form VI — A = (1 — y) has been intro- 
duced to remove leading singularities as A — ► 1. The re- 
sults for the spin- half model are shown in Fig. [T] and the 
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TABLE I: Series expansion coefficients for the spin-stiffness 
of the spin- 1/2 Heisenberg model on the anisotropic square 
lattice for selected values of the anisotropy coupling J y . 



order J x 


Jy 


Pan: 


Pay 


1 


0.9 


0.2500000000 


0.2250000000 


1 




0.0892857142 


0.0698275862 


2 




-0.0927175942 


-0.0926628428 


3 




-0.0151861349 


-0.0132149279 


4 




-0.0045297105 


-0.0011224771 


5 




0.0003357042 


0.0021422839 







-0.0057755038 


-0.0049661531 


7 




-0.0003145877 


-0.0002219010 


8 




-0.0042654841 


-0.0038964097 


1 


0.7 


0.2500000000 


0.1750000000 


1 




0.1041666666 


0.0453703703 


2 




-0.0868950718 


-0.0850367899 


3 




-0.0172969207 


-0.0105492432 


4 




-0.0045937341 


0.0050836143 


5 




-0.0012261711 


0.0039781766 


6 




-0.0059912113 


-0.0040000116 


7 




-0.0006353314 


-0.0005917251 


8 




-0.0045520522 


-0.0033419898 


1 


0.5 


0.2500000000 


0.1250000000 


1 




0.1250000000 


0.0250000000 


2 




-0.0891666666 


-0.0779166666 


3 




-0.0235044642 


-0.0081537698 


4 




0.0013212991 


0.0140706091 


5 




-0.0017270061 


0.0056040091 







-0.0060489454 


-0.0045071391 


7 




-0.0009351484 


-0.0017589857 


8 




-0.0057105087 


-0.0027490512 


n i 
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0.1562500000 


0.0097826086 


2 




-0.1094346417 


-0.0666305597 


3 




-0.0408511573 


-0.0052660964 


4 




0.0234241580 


0.0242642206 


5 




0.0016805213 


0.0056385893 


6 




-0.0085614425 


-0.0079166392 


7 




-0.0010470671 


-0.0031407541 


8 




-0.0085663482 


-0.0003268242 



results for spin-one model are shown in Fig. [2 In the ID 
limit, the spin-stiffness constant is know to be 0.25 J from 
exact calculations by Shastry and Sutherland^. This 
value is clearly larger than the square-lattice case where 
p s ~ 0.18J. Our results are more accurate away from 
the ID limit, but they clearly appear to approach the ID 
limit in a smooth and monotonic manner. For the spin- 
one case it is known that the Neel order disappears at 
an anisotropy ratio of approximately J y /J x = 0.0 1 13 ' 14 . 
The transition should be in the universality class of the 
3D Heisenberg model. The spin-stiffness should vanish 
at the transition with an exponent of v = 0.7—. The fit 
shows that the data agrees well with these expectations. 

For the spin-half case we also fit the small anisotropy 
regime to a linear behavior. The results arc shown in 
Fig. [3] where they are compared to the spin- wave results 
discussed in the next section. We find that the anisotropy 



TABLE II: Series expansion coefficients for the spin-stiffness 
of the spin-1 Heisenberg model on the anisotropic square lat- 



tice for selected values of the anisotropy couplin 


g Jy 


order J x J y 


Pax 


Pay 


1 0.9 


1.0000000000 


0.9000000000 


1 


0.1515151515 


0.1208955223 


2 


-0.1535289080 


-0.1422381534 


3 


0.0207519434 


0.0193249175 


4 


-0.0428180052 


-0.0391942098 


5 


0.0072984862 


0.0068989716 


(i 


-0.0210286555 


-0.0190202422 


7 


0.0043677522 


0.0041430202 


1 0.7 


1.0000000000 


0.7000000000 


1 


0.1724137931 


0.0803278688 


2 


-0.1553008729 


-0.1192091787 


■ ', 


0.0196401862 


0.0152963896 


4 


-0.0421528844 


-0.0313929868 


5 


0.0066352559 


0.0055840695 


6 


-0.0210149617 


-0.0149157761 


7 


0.0040511283 


0.0034008244 


1 0.5 


1.0000000000 


0.5000000000 


1 


0.2000000000 


0.0454545454 


2 


-0.1688941361 


-0.0982465564 


3 


0.0182177620 


0.0106543293 


4 


-0.0421971280 


-0.0242928541 


5 


0.0050652115 


0.0041087839 


6 


-0.0210666888 


-0.0106538747 


7 


0.0035338724 


0.0025923026 


1 0.3 


1.0000000000 


3000000000 


1 


0.2380952380 


0.0183673469 


2 


-0.2069165600 


-0.0746342335 


■ ', 


0.0173382982 


0.0055030901 


4 


-0.0479905685 


-0.0189077310 


5 


0.0002743242 


0.0023851969 


6 


-0.0213911544 


-0.0061871738 


7 


0.0021597090 


0.0016832884 


can be expressed as 



Psx/Psy = 1 + K{J x /J y - 1) , (5) 

where k = 1.8. This deviates significantly from the naive 
expectation k = 1. 



II. SPIN- WAVE CALCULATION 

In this Section we consider the self-consistent version 
of spin-wave theor y 15 ' 16 i 17 . To apply this approach we 
subdivide the lattice into sublattices A and B and use 
the Dyson-Maleev representation for spin operators on 
each sublattice, 

S+ = V2Sa, , S? = S - a\a u ieA (6) 
S7 = V2S(a\ - — atajai) 
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and 

St = V2Sb], S? = -S + b\b i , ieB (7) 
Sr = V2S(h - ^b\bA), 

where a\, a%, and b\, bi are the Bose operators. Introduc- 
ing the operators 

B, l = l a l ieA (8) 

and decoupling the four-boson terms in the Hamiltonian 
into all pos sible two-boson combinations, we derive (see 
Refs. [I5|ll6|) 

Wsswt = ^2 J m( B \ B i- B l + 8 B i) , ( 9 ) 

i,S 

where 6 = x,y correspond to the nearest neighbour sites 
in the x and y directions, 

75 = ~S+ (aih+s) (10) 

are the short-range order parameters and 

S = (SieA) = — (^ies) 

is the sublattice magnetization. Diagonalizing the Hamil- 
tonian (J9j) one finds the self-consistent equations at T = 

75 = S + ojr cos k s 

k k 

« = s+i /2-e|. (ID 

k 

where the antiferromagnetic spin- wave spectrum has the 
form 

£ k = ^/rg-r* , (12) 

with 

T k = 2 ( J x j x cos k x + Jy^j y cos k y ) (13) 

and To = Tk^o ( w ^ assume here that the ground state is 
antifcrromagnetically ordered, otherwise a bosonic chem- 
ical potential fx ^ should be introduced in the disper- 
sion to fulfill the condition S = 0, see Refs. Hill). 
The parameters 75 are simply related to the spin corre- 
lation function at the nearest-neighbor sites by 

ls = \(S t S l+s )\ 1/2 . (14) 

The spin- wave stiffhes along the x and y axes is expressed 
through these parameters as 

p s s = JsSSjs (15) 



The results obtained according to Eqs. (fTTj) and (fTS"]) are 
shown in Figs. [TJ [51 While for 5=1 the spin- wave 
results are close to those obtained by series expansions 
over the entire parameter range, for S = 1/2 at large 
anisotropy, the two techniques give results for p sy which 
are qualitatively different. This difference is most proba- 
bly due to topological excitations within the half-integer 
spin chains, which are not considered by spin-wave the- 
ory. For S = 1, at the transition to the Haldanc phase, 
we find the stiffness exponents to be v x = 0.50 and 
v y = 1.07, in qualitative agreement with series expan- 
sions. 

For S = 1/2, a quantitative discrepancy between spin- 
wave theory and series expansions is visible already at 
small anisotropics. While series expansion and spin- wave 
stiffness along the weaker exchange couplings axis (J y ) 
remain very close to each other, a discrepancy arises 
from the stiffness along the stronger coupling direction. 
We find the coefficient in the linear fit ([5]) k$w = 1-3, 
somewhat lower than what found by series expansion, as 
shown in Fig [3] It is possible that part of the difference 
is due to numerical inaccuracies or high order effects in 
1/5. However, with increasing anisotropy the difference 
is not just quantitative. It becomes qualitative and it 
implies the onset of new physics for the spin-half case 
associated with the Berry phase terms^. 

III. PINNING 

We base our considerations on the theory of under- 
doped cuprates suggested in Ref. 0. According to this 
theory, the ground state of an underdoped uniformly 
doped cuprate is a spin spiral, spontaneously directed 
along the (1,0) or the (0,1) direction. At sufficiently 
small doping, x < x c , the spiral has a static compo- 
nent while at x > x c it is fully dynamic. Here x is the 
concentration of holes in CuC>2 plane. According to this 
picture, the electronic liquid crystal observed in Ref. [HHI 
is the mostly dynamic spin spiral which may still have 
some small static component. For Sr doped La2CuC>4, 
the value of the critical concentration is x c « 0.11, and 
for YBa2Cu306.+j, is x c « 0.09. The absolute value of 
the wave vector of the spin spiral (static or dynamic) is 
given by 

Q=^-x. (16) 

Ps 

Here p s ss 0.18J is the spin-stiffness of the initial Heisen- 
berg models, J ~ 130meV is the antiferromagnetic ex- 
change parameter of the model, and g is the coupling 
constant for the interaction between mobile holes and 
spin waves. We set the spacing of the tetragonal lattice 
equal to unity, so the the wave vector Q is dimension- 
less. To fit the neutron scattering experimental data to 
the position of the incommensurate structure in Sr doped 
single layer La2Cu04, we need to set g w J, and to fit 
similar data for double layer YBa2Cu306.+i,, we need to 
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set g fa 0.7 J. It is not clear yet why the values of g 
for these compounds are slightly different, but for pur- 
poses of the present work this difference is not impor- 
tant. The coupling constant g was calculated within the 
extended t-J mode l 19 i 20 . The result is g = Zt where t 
is the nearest site hopping matrix clement and Z is the 
quasiholc residue. It is known that in cuprates t w 3J 
and Z ?s 0.3. Thus the calculated value of the coupling 
constant, g w J, agrees well with that found by fitting 
of experimental data. The ground state energy of the 
spin spiral state consists of two parts. The spin spiral 
with the wave vector Q gives rise to the gain —gQ in the 
kinetic energy of a single hole. On the other hand the 
spiral costs the spin elastic energy p s Q 2 /2. So the total 
balance is E = p s Q 2 /2 — xgQ, and minimization with 
respect to Q gives the wave vector (TTS)) and the energy 
per elementary cell 



E = psQ 1 12 - xgQ 



2p s 



(17) 



There are also quantum corrections to this energy, 
but they are small and hence not important for our 
purposes^. Note, that Eq. (fT?"]) is valid for both x < x c 
and x > x c , assuming that x is not large. 

Up to now we have disregarded the anisotropy assum- 
ing a perfect square lattice. To analyze anisotropy in the 
spiral direction wc have to replace 



t->t(l±e) 



(18) 



where e is due to the lattice deformation, so t a = t(l + e) 
and tb = t(l — e). The antifcrromagnetic exchange, J oc 
t 2 /U, also becomes anisotropic, J a = J(l + 2e), Jb = 
J(l — 2e). Hence the spin-stiffness is replaced by p s — > 
p s (l ± 2kc), where k w 1.8 has been calculated above. 
Now we can see how the lattice deformation influences 
the spiral energy (JT7J) . In the case when Q is directed 
along the a* = (1,0) we have to replace in (fT7|) t — > t a 
and p s — > p sa ] and in the case when Q is directed along 
the b* = (0, 1) we have to replace in (JTTJ) t — > tb and 
p s —* Psb- Note that the quasiparticle residue Z is a scalar 
property and therefore it is independent of direction of 
Q. Altogether, with account of the anisotropy, the energy 
(TTT)) is replaced by 



E 



g 2 x 2 
^p7 



[1 T 2(« - l)e] 



(19) 



The minus sign corresponds to Q directed along the a* = 
(1,0) axis, and the plus sign corresponds to Q directed 
alone the b* = (0, 1) axis. Interestingly, without the spin- 
quantum-fluctuations effect (i. c. if k = 1) the anisotropy 
in energy disappears. 

Since a* < b* it is most natural to assume that t a > tb, 
This means that e > 0. This point is supported by the 
LDA calculation performed in Ref. [2l|. In this case, ac- 
cording to Eq. (TT9"!) ■ the energy of the state with Q along 
the a*-axis is higher than that with Q along the 6*-axis. 



This disagrees with the experimental data in Ref. [l|[2|. 
However, the anisotropy of the hopping matrix element 
t is not straightforward. There are two competing con- 
tributions to e. The first one is related to the lattice 
deformation and is positive. The second one is related to 
oxygen chains that are present in YBa2Cu3 06.45 and is 
negative. In principle it is possible for the negative con- 
tribution to win, making e negative^. The neutron scat- 
tering anisotropy has been previously discussed within 
the Pomeranchuk instability scenario^. This is probably 
not sufficient to explain the newest data, see discussion 
in Ref. l24l. However, it is interesting to note that, to ex- 
plain the sign of the pinning, the Pomeranchuk scenario 
also requires a negative e. Anyway, for further numerical 
estimates, we will assume 



-0.02 



(20) 



The absolute value is consistent with the 1% lattice de- 
formation and the sign has been discussed above. In this 
case, according to Eq. (TT9")) . the energy of the state with 
Q along the a*-axis is lower and this is consistent with 
experimental data. The direction of pinning energy at 
x = 0.09 reads 



C(K-1) 



2 2 

9 x 



5 x W- 2 meV 



(21) 



This is the pinning energy per Cu site and it is a 
pretty strong pinning. For comparison, the pinning en- 
ergy of spin to the orthorhombic b-direction in undoped 
La2Cu04 is just ~ 1.5 x 10 _3 meV r . Assuming that the 
correlation length is at least comparable with the period 
of the spin spiral, £ ~ 2ir/Q ~ 17, we find that the total 
pinning energy per correlation unit is ~ 5 x 10~ 2 £ 2 ~ 
\bmeV , which is a significant energy scale. 



IV. CONCLUSIONS 

In this paper, we have studied the spin-stiffness con- 
stants for spatially anisotropic spin-half and spin-one 
Heiscnbcrg models using scries expansions and self- 
consistent spin-wave theory. The theoretical results have 
been of interest in themselves and show the importance 
of Berry phase interference terms in anisotropic square- 
lattice models. 

Our primary motivation for the study has been to un- 
derstand the phenomena of electronic liquid crystal and 
its pinning in high temperature superconductors. Wc 
find that quantum interference effects significantly en- 
hance the spin-stiffness anisotropy and this can provide 
the primary mechanism for the pinning of the liquid crys- 
tal direction. Wc have provided a detailed quantitative 
account of the pinning energy in YBCO. 



6 



Acknowledgments sions and comments. 



We are grateful to 0. K. Andersen, V. Hinkov, 
B. Keimer, and H. Yamasc for very important discus- 



1 V. Hinkov, P. Bourges, S. Pailhes, Y. Sidis, A. Ivanov, C. 
D. Frost, T. G. Perring, C. T. Lin, D. P. Chen, B. Keimer, 
Nature Physics 3, 780 (2007). 

2 V. Hinkov, D. Haug, B. Fauque, P. Bourges, Y. Sidis, 
A. Ivanov, C. Bernhard, C. T. Lin, and B. Keimer, Sci- 
ence 319, 597 (2008). 

3 I. Affleck and B. I. Halperin, Journal of Physics A: Math- 
ematical and General 29, 2627 (1996). 

4 I. Affeck, M. P. Gelfand, and R. R. P. Singh, Journal of 
Physics A: Mathematical and General 27, 7313 (1994). 

5 A. Parola, S. Sorella, and Q. F. Zhong, Phys. Rev. Lett. 
71, 4393 (1993). 

6 A. W. Sandvik, Phys. Rev. Lett. 83, 3069 (1999). 

7 V. Y. Irkhin and A. A. Katanin, Phys. Rev. B 61, 6757 
(2000). 

8 B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 65, 243 
(1990). 

9 S. Sachdev, Low Dimensional Quantum Field Theories for 
Condensed Matter Physicists (World Scientific, Singapore, 
1995), edited by Y. Lu, S. Lundqvist, and G. Morandi. 

10 R. R. P. Singh and D. A. Huse, Phys. Rev. B 40, 7247 
(1989). 

11 C. J. Hamer, Z. Weihong, and J. Oitmaa, Phys. Rev. B 



50, 6877 (1994). 

12 J. Oitmaa, C. Hamer, and W. Zheng, Series Expansion 
Methods for Strongly Interacting Lattice Models (Cam- 
bridge University Press, New York, 2006). 

13 I. Affleck, Phys. Rev. Lett. 62, 474 (1989). 

14 T. Pardini and R. R. P. Singh, arXiv:0801.3855vl [cond- 
mat.str-el]. 

15 V. Y. Irkhin, A. A. Katanin, and M. I. Katsnelson, Phys. 
Rev. B 60, 1082 (1999). 

16 A. A. Katanin and V. Y. Irkhin, Physics-Uspekhii 50, 613 
(2007). 

17 M. Takahashi, Phys. Rev. B 40, 2494 (1989). 

18 A. I. Milstein and O. P. Sushkov, arXiv:0711.4865. 

19 J.-i. Igarashi and P. Fulde, Phys. Rev. B 45, 10419 (1992). 

20 O. P. Sushkov and V. N. Kotov, Phys. Rev. B 70, 024503 
(2004). 

21 O. K. Andersen, A. I. Lichtenstein, O. Jepsen, and 
F. Paulsen, J. Phys. Chem. Solids 56, 1573 (1995). 

22 O. K. Andersen, private communication. 

23 H. Yamase and W. Metzner, Physical Review B 73, 214517 
(2006). 

24 H. Yamase, arXiv:0802.1149 [cond-mat.str-el]. 



